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Abstract 

Using first-principles density-functional calculations, we explore the possibility of magnetic order at the 
rebonded D# step of the Si(001) surface. The rebonded Db step containing threefold coordinated Si atoms 
can be treated as a one-dimensional dangling-bond (DB) wire along the step edge. We find that Si atoms 
composing the step edge are displaced up and down alternatively due to Jahn-Teller-like distortion, but, if 
Si dimers on the terrace are passivated by H atoms, the antiferromagnetic (AFM) order can be stabilized 
at the step edge with a suppression of Jahn-Teller-like distortion. We also find that the energy preference 
of AFM order over Jahn-Teller-like distortion is enhanced in an oscillatory way as the length of DB wires 
decreases, showing the so-called quantum size effects. 
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I. INTRODUCTION 



Creation of a magnetically ordered state in carbon materials has attracted much attention be- 
cause of its fundamental scientific interest as well as for the potential applications in diverse fields.- 
Indeed, a number of theoretical and experimental studies have been performed to explore the pos- 
sibility of magnetic orders in various carbon materials such as graphite,- carbon nanofoams,- 
nanodiamonds, 4 fullerenes,- 1 ^ and graphene. 7- - For instance, it was theoretically predicted that 
graphene flakes or nanoribbons with zigzag edges have a high density of states at the Fermi level, 
therefore yielding a ferromagnetic order at the edge.-"- 

Since silicon belongs to the same group IV elements with carbon, it is very interesting to extend 
the discussion of magnetic order in silicon materials. Recently, a first-principles density-functional 
theory (DFT) calculation of Erwin and Himpsel 10 predicted that Au-induced vicinal Si surfaces 
such as Si(553)-Au and Si(557)-Au exhibit an antiferromagnetic (AFM) coupling between the 
magnetic moments arising from the Si dangling-bond (DB) electrons. Here, the critical structural 
motif for the presence of AFM order is a graphitic honeycomb ribbon structure of Si atoms at the 
step edge. In this sense, the novel way of manipulating magnetism in silicon materials is to create 
one-dimensional (ID) nanostructures at step edges, which provide more favorable conditions for 
long-range structural order with an array of spins. 

In this paper, we explore the possibility of magnetic orders at the stepped Si(OOl) surface using 
first-principles DFT calculations. Here, we consider the double-atomic height step, i.e., the so- 
called rebonded Dg step [see Fig. 1(a)], which prevails on the Si(OOl) surface.— This step contains 
threefold coordinated Si atoms, giving rise to a ID array of DBs along the step. Our spin-polarized 
DFT calculations for such a ID DB wire show that the AFM or ferromagnetic (FM) configuration, 
where neighboring spins of DB electrons are coupled antiferromagnetically or ferromagnetically 
with each other, is energetically unfavored over the nonmagnetic (NM) configuration, where Si 
atoms composing the rebonded Dg step are alternatively displaced up and down due to Jahn- 
Teller-like distortion. However, we find that, if Si dimers on the terrace are passivated by H 
atoms, DB electrons at the step become more localized to yield the preference of the AFM order 
with a suppression of Jahn-Teller-like distortion. In addition, we find that, as the length of DB 
wires decreases, the energy preference of the AF configuration over the NM one increases in an 
oscillatory way. Our findings demonstrate that the AFM spin ordering of DB electrons at the 
rebonded Dg step can be stabilized by passivating the terrace with atomic H, and its strength can 
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be manipulated with respect to the length of DB wires. 




FIG. 1: (color online) Optimized structures for the rebonded Db step of a Si( 1 17) surface within (a) the NM 
configuration and (b) the AFM configuration. The large and small circles represent Si and H atoms, respec- 
tively. For distinction, Si atoms at the step edge (denoted by S) and Si dimers on the terrace (denoted by Ti 
and T2) are drawn in different dark colors. The optimized structures of an infinitely long DB wire within 
the NM and AFM configurations as well as a finitely long DB wire (DB-2) within the AFM configuration 
are also given in (c), (d), and (e), respectively. The insets in (a) and (e) show the top view and the side view 
of each structure, respectively. 



II. CALCULATION METHOD 

The total-energy and force calculations were performed by using nonspin-polarized and spin- 
polarized DFT— ^ within the generalized-gradient approximation.— The Si and H atoms were 
described by norm-conserving pseudopotentials.— To simulate the rebonded Dg step, we used a 
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triclinic supercell with the lattice vectors of a = (\/2, 0, 0), b = (-1/2 a/2, 1/2^2, 1/2), and c = (0, 
0, 17/4) in units of the optimized Si lattice constant ao = 5.480 A [see Fig. 1(a)], corresponding to 
a vicinal surface with Miller indices (117). The stepped Si(001) surface was modeled by a peri- 
odic slab geometry containing a thickness of six Si atomic layers along the [001] direction. In this 
periodic slab geometry, periodically repeated rebonded Dg steps were separated by terraces with a 
width of two Si dimers [see Fig. 1(a)]. The bottom Si layer of the slab was fully passivated by one 
or two H atoms per Si atom. For simulations of finitely long DB wires, we used a supercell with 
the lattice vectors of 4a, b, and c, which includes eight Si atoms along the step edge. The electronic 
wave functions were expanded in a plane-wave basis set with a cutoff of 20 Ry, and the k-space 
integration was done with 4x2x1 (1x2x1 for finitely long DB wires) k-points in the surface Bril- 
louin zone.— All atoms except the bottom two Si layers were allowed to relax along the calculated 
Hellmann-Feynman forces until all residual force components were less than 1 mRy/bohr. The 
present calculation scheme has been successfully applied not only for the electronic and magnetic 
properties of various DB wires on the H-passivated C(001), Si(001), and Ge(001) surfaces^ but 
also for the reaction and adsorption of various hydrocarbons on the Si(100) surface.— 

III. RESULTS 

We begin to optimize the atomic structure of the rebonded Dg step of the Si(001) surface. 
Figure 1(a) shows the optimized structure of the rebonded Dg step, where the step edge atoms 
are designated as S. The rebonded D B step separates upper and lower terraces containing two 
buckled dimers T\ and T2 [see Fig. 1(a)]. We find that S atoms have relatively long bonds to 
the step atoms (located between S and T\) with two different bond lengths of 2.46 and 2.55 A. 
Furthermore, S atoms are displaced up and down alternatively with a height difference of 0.66 
A, in good agreement with that (0.67 A) obtained from a previous 1 - DFT calculation. The dimer 
bond length Jdimer and the dimer height difference Ah of T\ are found to be 2.36 and 0.71 A, 
respectively, significantly different from those (2.29 and 0.49 A) of Tj. The relatively smaller 
value of Ah in T2 compared with T\ is possibly due to a strain effect at the step edge: i.e., half of 
the subsurface atoms attached to Ti retain bulk-like bonding arrangements so that T2 has too much 
strain for buckling. We note, however, that Jdimer and Ah of T\ are close to those (2.36 and 0.71 
A) of buckled dimers at the flat Si(001) surface. 

Each S atom at the Dg step has one unpaired electron, forming a ID array of DBs along the 
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step edge. This kind of DB wire is similar to the case fabricated by the selective removal of H 
atoms from an H-passivated Si(OOl) surface along one side of an Si dimer rowM^- For the latter 
DB wires of finite lengths, recent DFT calculation s 17 ' 23 predicted that the AFM configuration is 
energetically favored over the NM one and the energy difference between the two configurations 
varies with respect to the wire length. 24 In order to examine the possibility of magnetic orders at 
the Dg step, we determine the atomic structure of the Dg step within the AFM configuration. The 
optimized structure of the AFM configuration is shown in Fig. 1(b). Unlike the NM configuration 
showing a Jahn-Teller-like distortion, the AFM configuration is found to suppress the up and 
down displacements of S atoms, indicating the absence of spin-Peierls instability. We find that 
the NM configuration is more stable than the AFM configuration by 28 meV/DB. This result is 
different from our previous^ 2 - result for the DB wire formed on the H-passivated Si(OOl) surface, 
where the AFM configuration is slightly favored over the NM configuration by 8 meV/DB. This 
different energetics of the NM and AFM configurations between the Dg step and the H-passivated 
Si(OOl) surface may be attributed to more delocalized character of DB electrons at the Dg step, as 
discussed below. 

Figure 2(a) shows the electronic band structure of the Dg step, together with the charge char- 
acters of three surface states (55i, SS2, and SS3) located below the Fermi level. It is seen that the 
three surface states have different charge characters depending on the k vectors. For instance, the 
charge character of 55 1 at the T (X) point represents the DB state of T2 (T\), indicating a hybridiza- 
tion of the DB states originating from T\ and T2. In this way, 552 (553) involves a hybridization of 
the DB states originating from 5 and T\ (5 and T2). Thus, we can say that DB electrons at the Dg 
step are much delocalized to favor the NM configuration rather than the AFM configuration. 

Recently, several theoretical studies predicted surface magnetism induced by the adsorption 
of hydrogen atoms on surfaces.— Amanna et a/.— reported an H-induced magnetic order 
on the Pd(llO) surface. Okada et al.— predicted a ferrimagnetic spin ordering in various DB 
networks created on an H-passivated Si(lll) surface. Similarly, Lee et al.— predicted an AFM 
spin ordering in DB wires fabricated on an H-passivated Si(OOl) surface. On the basis of these 
previous studies, we speculate that the derealization of DB electrons at the Dg step could be 
prevented by passivating the terrace with atomic H [see Fig. 1(c) and 1(d)], thereby giving rise 
to the stabilization of magnetic order: that is, the H-passivation of terrace gives a confinement 
of DB electrons at the step edge, producing magnetic order due to an enhanced electron-electron 
interaction. The fabrication of such a DB wire can be realized by passivating a stepped Si(OOl) 
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surface with atomic H and then removing H atoms along the step edge via the scanning tunneling 
microscope nanolithography technique.—"— This technique has been utilized to fabricate the DB 
wires on the H-passivated Si(OOl) surface.—^ 

We optimize the atomic structure of an infinitely long DB wire (designated as DB-°o) formed 
on the Dg step within the NM, FM, and AFM configurations. We find that the NM configura- 
tion shows alternating up and down displacements of S atoms with a height difference of 0.77 A, 
whereas the FM and AFM configurations have a nearly unbuckled structure. The AFM configu- 
ration is found to be more stable than the NM (FM) configuration by 14 (48) meV/DB. Here, the 
total spin of the AFM configuration is zero, while that of the FM configuration is 1 /j#/DB. Figure 
2(b), 2(c), and 2(d) show the calculated band structures of the NM, FM, and AFM configurations, 
respectively. The charge characters of the DB state (designated as SS) reveal that the NM con- 




FIG. 2: (color online) (a) Calculated band structure of the rebonded Db step, together with the charge 
characters of three surface states (SS\ , SS2, and SS3) at the T and X points. The first contour line and the line 
spacings are the same as 0.001 electron/bohr 3 . The energy zero represents the Fermi level. The direction of 
F-X line is along the step. The results for an infinitely long DB wire created by the H-passivation of terrace 
are also given in (b), (c), and (d) for the NM, FM, and AFM configurations, respectively. In (c) and (d), the 
solid (dashed) lines represent the subbands of up (down) spin. Charge contour plots in (b), (c), and (d) are 
drawn in a vertical plane containing S atoms. 
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figuration has a charge transfer from the down to the up S atoms [see Fig. 2(b)] while the AFM 
configuration involves an opposite spin orientation between adjacent DBs [see Fig. 2(d)]. We note 
that the band gap of the AFM configuration is 0.64 eV, larger than those (0.31 and 0.38 eV) of the 
FM and NM configurations. This manifests that the electronic energy gain by AFM spin ordering 
is greater than that by FM spin ordering or Peierls-like lattice distortion. 

It was reported that in the DB wire fabricated on an H-passivated Si(001) surface, electron- 
electron correlations are competing with electron-lattice couplings due to Peierls instability or 
Jahn-Teller effect.— To find the on-site electron-electron interactions in the DB wire [Fig. 
1(d)] formed on the Dg step, we evaluate the Hubbard correlation energy U by using the con- 
strained DFT calculations. 30 Here, we simulate the transfer of one spin-down electron from the 
SS state [in Fig. 2(d)] to the lowest unoccupied state and calculate the change in the total energy. 
We obtain U = 1.71 eV, which is much larger than the electron hopping parameter t = 0.23 eV 
(estimated from the bandwidth of the SS state). The resulting large ratio of U/t = 1.43 for an 
infinitely long DB wire indicates the importance of electron correlations, giving rise to the prefer- 
ence of an AFM spin ordering over a Peierls-like distortion. From the energy difference between 
the FM and AFM configurations, the exchange coupling constant / between adjacent spins can be 
estimated as 48 meV. Note that S atoms with highly localized DB electrons are covalently bonded 
to the subsurface Si atoms which have a negligible spin polarization. Thus, it is most likely that 
an individual DB is fully spin-polarized by an intra-atomic exchange, and such magnetic moments 
of S atoms are antiferromangnetically coupled with each other via the superexchange mechanism, 
similar to the case of the DB wires formed on the H-passivated Si(001) surface.— 22£L 

In order to see how the magnetic stability of the DB wires varies with respect to the wire 
length, we determine the atomic structures of finitely long DB wires within the NM, FM, and AFM 
configurations. Here, we consider finitely long DB wires containing up to six DBs (designated as 
DB-2, DB-3, DB-4, DB-5, and DB-6, respectively). The calculated total energies of the NM and 
FM configurations relative to the AFM configuration as a function of the wire length are plotted 
in Fig. 3. We find that the AFM configuration is more stable than the NM configuration by 
A£ NM -afm = 89, 108, 68, 81, and 51 meV for DB-2, DB-3, DB-4, DB-5, and DB-6, respectively. 
Thus, the energy preference of AFM ordering over Jahn-Teller-like distortion tends to increase 
as the wire length decreases (see Fig. 3). This may be ascribed to the fact that the two ends in 



7 



120 




20 - \ 



n I i i i i i i 

2 3 4 5 6 CO 

DB-wire length 

FIG. 3: (color online) Calculated total energies of the NM and FM configurations relative to the AFM 
configuration for finitely long DB wires as a function of the wire length. The results for DB-oo are also 
given. 

shorter DB wires give rise to a relatively larger elastic-energy cost for the alternating up and down 
displacements in the NM configuratin. It is noteworthy that A£nm-afm for DB3 (DB5) is larger 
than those of DB2 and DB4 (DB4 and DB6). Such an even-odd oscillatory behavior correlates 
with the variation of DB states with respect to the wire length, as discussed below. Figure 3 shows 
that the AFM configuration is energetically favored over the FM one by 29, 33, 32, 38, 43, and 
48 meV/DB for DB-2, DB-3, DB-4, DB-5, DB-6, and DB-°o, respectively. Using the Heisenberg 
spin Hamiltonian H = J^ij^i ' Sj where we assume only the existence of the nearest neighbor 
interactions in DB wires, we obtain J = 57, 49, 42, 47, 52, and 48 meV for DB-2, DB-3, DB-4, 
DB-5, DB-6, and DB-oo, respectively. 31 This result indicates that the exchange coupling constant 
between adjacent spins varies slightly with respect to the wire length. 



TABLE I: Calculated band gaps (in eV) for the DB-2, DB-3, DB-4, DB-5, DB-6, and DB-°o wires within 
the NM and AFM configurations. 





DB-2 


DB-3 


DB-4 


DB-5 


DB-6 


DB-oo 


NM 


0.43 


no 


0.23 


no 


0.24 


0.38 


AFM 


0.94 


0.89 


0.87 


0.82 


0.80 


0.63 



Figure 4(a), 4(b), 4(c), and 4(d) show the band structures for DB-2, DB-3, DB-4, and DB-5, 
respectively. It is seen that the subbands due to DB electrons exist near Ep . We find that, as the 
wire length increases, the subbands are lower in energy and their spacings decrease, showing the 
evolution of quantum well states. For DB-3 and DB-5, the NM configuration has a half-filled 
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FIG. 4: (color online) Calculated band structures of (a) DB-2, (b) DB-3, (c) DB-4, and (d) DB-5 within 
the NM and AFM configurations. The energy zero represents the Fermi level. The direction of F-X line is 
along the DB wire. The lines represent the subbands due to DB electrons. In the AFM configuration, the 
subbands of up (down) spin are drawn with the solid (dashed) lines. 

subband, while the AFM configuration has a band-gap opening of 0.89 and 0.82 eV, respectively 
(see Table I). On the other hand, for the even-numbered DB wires of DB-2, DB-4, and DB-6, 
the NM configuration has a band gap of 0.43, 0.23, and 0.24 eV, respectively, which are smaller 
than the corresponding ones (0.94, 0.87, and 0.80 eV) of the AFM configuration. These different 
features of subbands between odd- and even-numbered DB wires represent the so-called quantum- 
size effects, accounting for the above-mentioned oscillatory behavior of A£nm-af with respect to 
the wire length (see Fig. 3). The calculated band gaps for finitely long DB wires also manifest that 
the electronic energy gain by the AFM spin polarization is larger than that by the Jahn-Teller-like 
lattice distortion. 

IV. SUMMARY 



Our spin-polarized DFT calculations predicted that the rebonded Dg step of the Si(001) sur- 
face can stabilize an AFM order by passivating the terrace with H atoms. We found that the 
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H-passivation of terrace enhances the localization of DB electrons at the step edge, leading to an 
AFM coupling between adjacent spins of DB electrons. We also found that finitely long DB wires 
exhibit the evolution of quantum-well states as a function of the wire length, giving rise to the 
even-odd oscillatory behaviors in their band structures and energetics. The H-induced magnetism 
predicted here in the Si DB wires may offer a unique way to device the spintronics on semiconduc- 
tor substrates. Finally, we anticipate that scanning tunneling microscopy can identify indirectly 
the predicted AFM order by measuring the height profile along the DB wire, because the AFM 
configuration involves no buckling of the Si atoms composing the DB wire. 
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